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We show that the number of observed voids in galaxy redshift surveys is a sensitive function of 
the equation of state of dark energy. Using the Fisher matrix formalism we hnd the error ellipses 
in the wo — Wa plane when the equation of state of dark energy is assumed to be of the form 
wcpl[z) = uiq + Waz/il + z). We forecast the number of voids to be observed with the ESA Euclid 
satellite and the NASA WFIRST mission, taking into account updated details of the surveys to 
reach accurate estimates of their power. The theoretical model for the forecast of the number of 
voids is based on matches between abundances in simulations and the analytical prediction. To 
take into account the uncertainties within the model, we marginalize over its free parameters when 
calculating the Fisher matrices. The addition of the void abundance constraints to the data from 
Planck, HST and supernova survey data noticeably tighten the wo — Wa parameter space. We thus 
quantify the improvement in the constraints due to the use of voids and demonstrate that the void 
abundance is a sensitive new probe for the dark energy equation of state. 

PACS numbers: 


I. INTRODUCTION 

The measurement of luminosity distances from super¬ 
novae of type la [HE] gives observational evidence for 
the accelerated expansion of the Universe. To account 
for such acceleration, the standard cosmological model 
relies on the existence of a mysterious component: dark 
energy. Different models of dark energy are viable, such 
as models where the Universe is permeated by a con¬ 
stant energy density (e.g. a cosmological constant with 
wq = —1) or models of dynamical dark energy. A pre¬ 
cise measurement of the equation of state of dark energy 
allows to distinguish among such models. 

The large scale structure of the Universe depends on 
the cosmological model, thus its study allows to trace 
the rate of expansion of our Universe. For instance, the 
number of galaxy clusters depends on the properties of 
dark energy: a comparison of observational counts with 
analytical predictions permits to constrain the dark en¬ 
ergy equation of state |3]. While this measure relies on 
the high density regions of the universe, the use of the 
low density ones, such as cosmic voids, is also sensitive to 
the properties of dark energy. Additionally, while clus¬ 
ter abundances can be problematic because of the high 
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intrinsic scatter in their mass-observable relationships, 
voids promise to suffer from this uncertainty in a milder 
(and different) way, as their volumes can be observed and 
compared to theory directly. 

Cosmic voids have remained fairly unexplored until re¬ 
cent times, due to the difficulty in acquiring sufficient 
data from sparsely populated regions of the Universe. 
Nowadays, as modern surveys such as the SDSS map out 
the cosmic web in great detail, entire void catalogs are 
being compiled HI] and used in a wide range of cosmo¬ 
logical applications, spanning from the Alcock-Paczyhski 
test [HIE] to the measurement of the integrated Sachs- 
Wolfe effect |3]. Analogously to galaxy cluster number 
counts, here we propose to use void number counts to 
constrain the cosmological model: the number of ob¬ 
served voids is sensitive to the dark energy equation of 
state. This means that simply counting voids in galaxy 
surveys enables us to constrain the properties of the lat¬ 
ter. 

Theoretical models based on the Press-Schechter for¬ 
malism m adapted for the case of voids (known as the 
Sheth and Van de Weygaert model m) and the recent 
observational work with voids allows in principle to fore¬ 
cast the number of voids that should be observed by a 
particular survey and with a given cosmological model. 
While recent works have extended the original Sheth and 
Van de Weygaert model (such as [EHIl]), for this work 
we will consider the original model, in the framework 
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of which we are qualitatively able to match simulations 
mimicking the tracer number density of upcoming sur¬ 
veys (as explained throughout the paper). Thus, the 
predicted abundance of observed voids for each galaxy 
redshift survey gives the power of the survey to constrain 
the dark energy equation of state through void counting. 
To reach this goal, we use the Fisher matrix formalism 
[mils] (as it has been done for different probes by m 
and [IH]). 

In this paper we focus on two upcoming surveys, Euclid 
[19] and WFIRST (Wide Field Infrared Space Telescope) 
[iniisi], to establish their expected power to constrain 
dark energy in terms of void abundance. The two sur¬ 
veys are in some sense complementary. The Euclid satel¬ 
lite is a Medium Class mission of the ESA-led Cosmic 
Vision 2015-2025 that will cover 15000 deg^ of the sky. 
It will give pseudo-spectroscopic redshifts of tens of mil¬ 
lions of galaxies up to high redshift. On the other hand, 
WFIRST — a NASA-led mission — will cover a smaller 
portion of the sky, 2000 deg^, but with a higher sam¬ 
pling density of galaxies. The two surveys thus present 
two complementary strategies, on the one side a wide but 
shallower survey (Euclid); on the other side a deeper but 
narrower field (WFIRST). 

Naturally, one might ask which approach can provide 
the largest number of voids. Because cosmic voids exhibit 
a hierarchical structure HUES], an increase in the density 
of galaxies (a higher resolution sampling of the large- 
scale structure of the Universe) allows resolving voids of 
smaller size at each redshift. On the other hand, a wider 
field also increases the number of observed voids, so it 
is not obvious which survey strategy is best suited for 
maximizing the number of observed voids. 

In this paper we aim to answer this question by pre¬ 
senting a forecast for the number of voids to be detected 
by Euclid and WFIRST. We also establish the constrain¬ 
ing power from their void abundances in view of the 
dark energy equation of state. To this end, we consider 
the Chevallier-Polarski-Linder [23l [24] parametrization 
for the dark energy equation of state: 

2 ) 

WCPLiz) = Wq +Wa - — (1) 

Z+1 

Once the number of observed voids is known, the 
Fisher matrix formalism can provide error ellipses in the 
WQ — Wa plane. The Fisher formalism gives an estimate of 
the constraints to be expected from cosmological probes 
in terms of area in parameter space. 

Recent work on observational data (e.g. the SDSS DR7 
and DR9 data, see I3ES]) sheds more light on the be¬ 
haviour of voids that are defined in the distribution of 
galaxies (as opposed to dark matter [2S1 EZj), thereby 
allowing a more robust estimation of void abundances 
based on the use of simulations and mock galaxy cata¬ 
logs that are calibrated to reproduce the properties of the 
data in the most realistic way possible. 

A number of cosmological tests based on void statis¬ 
tics have been proposed, e.g. Alcock-Paczyhski test on 


void stacks [TjHIHlIMj, tests of coupled dark energy- 
dark matter models [29] , tests of modified gravity models 
[5nH5^ , constraints from CMB lensing with cosmic voids 
[55] . To understand the potential of void-based cosmol¬ 
ogy with Euclid and WFIRST requires realistic estimates 
of void abundances based on updated details on both mis¬ 
sions, such as those we present in this paper. 

The paper is organized as follows: in Section 2 we 
present the theoretical calculation of the void abundance. 
In Section 3 we discuss the determination of the mini¬ 
mum void radius to be considered useful for our analysis. 
Section 4 presents the Fisher matrix analysis and in Sec¬ 
tion 5 we describe the simulations used for this work. 
Finally in Section 6 we present our results and conclude 
in Section 7. 


II. VOID NUMBER CALCULATION 

In order to estimate the constraining power from voids 
that can be observed with future surveys, we need to 
have a theoretical prediction of their expected abun¬ 
dance. This, in turn, requires a theoretical model for 
the formation of voids. One possibility to define the for¬ 
mation of a void is to consider the spherical expansion of 
a top-hat perturbation, analogously to the case of halo 
formation by spherical collapse [M], but with opposite 
sign. The moment of void formation can be defined as 
the evolutionary stage where two adjacent radial shells 
start crossing each other, leading to the formation of a 
ridge m- This consideration allowed Sheth and Van de 
Weygaert m to construct a model for the distribution 
of void sizes. We take this model as a starting point for 
determining void abundances. 

Following the Press-Schechter formalism m, we cal¬ 
culate the abundance of voids using 

M^n(M.») dM ^ ^ 

p M ' V ^ ' 

where M is the void mass, p the background density, and 
n{M, z) the number density of voids of given mass and 
redshift. According to Sheth and Van de Weygaert 
assuming Gaussian initial conditions, the fraction /(zz) 
of mass that has evolved into voids can be approximated 
as: 

where v = is the critical underdensity of void 

formation, linearly extrapolated to the time when shell¬ 
crossing takes place at the void edge, and cr^(M) is the 
variance of linear density fluctuations filtered on a scale 
R = (ff|)~^'^^- We anticipate here that we consider 
equation [^ in the case were the two-barrier distribution 
reduces to a single barrier, a particular case described 
in m where the void-in-cloud process is neglected, as 
detailed in the next Section. 
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Recent work measured the abundance of voids in 
galaxy surveys [e.g. m) and compared it to mock cata¬ 
logs from simulations with properties matching the data. 
The Sheth and van de Weygaert excursion set model can 
be used to model the void abundances accurately, as long 
as the value of Sy is taken as a free parameter. While 
this is a phenomenological approach, it uses the power 
of the theoretical description combined with the added 
realism by fitting a parameter to simulations matching 
observations. Indeed the distribution of voids needs to 
be validated with simulations in order to be reliable m- 

Thus, in the redshift range of interest to observe voids, 
we have used sub-sampled simulations matched to the 
tracer density of Euclid and WFIRST to tune this pa¬ 
rameter optimally. The value of Sy obtained with this 
methodology is —0.45, we find this value to be stable 
when varying the cosmological parameters. The value of 
Sy differs from the Sheth and van de Weygaert prediction 
because watershed voids are not spherical and do not ex¬ 
hibit uniform shell-crossing. As mentioned, recent works 
show that some improvements can be obtained |13) by 
validating the void distribution with mocks. The use of 
simulations tuned to observed data [5^ thus enhances 
the robustness of the abundance estimation we perform. 

Additionally, since the value of Sy might partially de¬ 
pend on the void definition, we use the same void finder 
(vide [3S], based on ZOBOV [3S]) for the treatment of ex¬ 
isting data and simulations with similar properties — to 
minimize bias due to the void definition. This is only 
a consistency choice, since generally, we do not expect 
void abundances to change much for void finders based 
on a similar concept as VIDE (namely a tessellation of the 
tracer distribution followed by a watershed transform to 
detect the structure of the cosmic web). 

In Equation!^ we used a, the linearly extrapolated vari¬ 
ance of density fluctuations S smoothed on the filtering 
scale i?Lag by the Fourier transform of the window func¬ 
tion for a top-hat filter W{x) = 3/x^ [sin(x) — a;cos(x)]. 
Considering the matter power spectrum Ps{k,a), we 
write the definition of a: 



k^Ps{k, a) 


WikRi^^giM)) 


2 M 
~k 


( 4 ) 


The calculation of <j depends on cosmology, thus it can 
be used to constrain cosmological parameters through the 
observed abundance of voids in a survey. We can change 
variables from the mass of the void to its radius (either 
Lagrangian or Eulerian) by using mass conservation in¬ 
side the void: 


— “^-^Lag P — (1 + ^v)-^-R euI P) (5) 

where Ay is the underdensity of the matter inside the 
void. From the spherical collapse model one finds Ay ~ 
—0.8 at the time of shell crossing. Since the expansion 
of the void reduces considerably after shell-crossing (see 


[TT1I33) we will assume this value of Ay for all voids irre¬ 
spective of when the shell formed around them. It should 
be noted that we have written the Eulerian radius in the 
comoving form, hence the background density is constant 
and equal to its present day value. To prevent confusion 
we will keep the Eulerian and Lagrangian subscripts for 
the radii throughout the paper. The observable quantity 
we measure in a galaxy redshift survey is the number of 
voids larger than a given radius located in a speci¬ 
fied redshift interval. It can be found as 


(•z+Az 


Ny_ = 


dz 


di^Eul 


dnn(i?Eui, z) 


dE 

dzdn 


( 6 ) 

where dV is the comoving volume element and the an¬ 
gular integration is over the angular size of the survey. 
The differential volume element is a function of cosmol¬ 
ogy through the Hubble rate. 

Considering Equations and the number of voids 
will depend on the cosmological model through a and the 
differential volume element. Additionally, the minimum 
radius of voids considered is a relevant quantity that will 
affect the forecast of void abundance. Its choice will be 
discussed in the next Section. 


III. MINIMUM VOID RADIUS 
DETERMINATION 

The constraining power of a void survey crucially de¬ 
pends on the minimum observable void radius. Particular 
care must be taken to determine this minimum radius for 
the voids in the survey: the abundance of voids increases 
rapidly with decreasing radius (analogously to the abun¬ 
dance of clusters). We consider two different criteria for 
defining the minimum void radius and take into account 
the most stringent [m, *.e. the largest radius of the 
two. 

The first and simplest criterion is based on the mean 
particle separation of galaxies in the survey. Taking into 
account the features of the survey, we consider as a mini¬ 
mum radius of a void twice the mean particle separation. 
Empty spaces among galaxies might not correspond to 
true voids in the dark matter distribution, as they may 
simply result from the discrete sampling of the underly¬ 
ing distribution (shot noise). 

Voids can be defined reliably down to the mean par¬ 
ticle separation, but considering twice the mean parti¬ 
cle separation 2Rmpa also guarantees a small impact of 
peculiar velocities in redshift space [3H]. Indeed, voids 
below that limit can be washed out or created by the 
effect of peculiar velocities. Therefore, to increase our 
confidence that an observed underdensity in the galaxy 
distribution is a true void and to limit systematic effects 
due to velocities, we demand that the radius of the void 
must be larger than twice the mean particle separation 
Rf°^ = 2.0 X (I/hgai)^/^ = 2.0 X 
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The comoving number density of galaxies figai depends 
on survey specifications, such as the maximum apparent 
luminosity observable or the galaxy luminosity function 
(see [T7] for more details). While the Euclid redshift sur¬ 
vey has undergone some changes recently (which might 
lead to updated values); to refer as much as possible to 
realistic values for the galaxy number densities of the sur¬ 
vey, for the present paper we will use as representative 
for each redshift bin the values reported by the Review of 
the Euclid Theory Working Group [39] . For the WFIRST 
satellite, while our results do not change much using 
the values in the Wide-field infrared survey telescope- 
astrophysics focused telescope assets WFIRST-AFTA fi¬ 
nal report [21j . we use updated values from the most re¬ 
cent Wide-Field InfrarRed Survey Telescope-Astrophysics 
Focused Telescope Assets WFIRST-AFTA 2015 Report 

m- 

A second threshold for the minimum observable void 
size is set by the void-in-cloud effect m- It is due to 
the fact that small voids tend to be located inside sur¬ 
rounding overdensities, and therefore will disappear with 
the gravitational collapse of the overdense region around 
them. In such cases the Press and Schechter formalism 
must be extended. As proposed by Sheth and Van de 
Weygaert m, a two-barrier excursion set model should 
be used — the need for two parameters for the void size 
distribution is unique to the void case. It involves solving 
a diffusion equation with two absorbing barriers instead 
of one, one for halos and one for voids. 

The void-in-cloud effect leads to a turn-over radius be¬ 
low which the number density of voids decreases towards 
smaller radii. For our analysis, as discussed in Section 2, 
we want to avoid this regime, so we consider voids with 
radii larger than the Lagrangian void-in-cloud radius for 
which we have , z) ~ 1. Such voids are unlikely 

to be embedded in an overdense region, hence giving a 
framework where the Press and Schechter formalism can 
be applied (only the smallest voids are subject to col¬ 
lapse m , [41] ). The analysis of velocity profiles confirms 
this scenario [42H44] . We thus calculate this radius, con¬ 
vert it from Lagrangian to Eulerian radius with equation 
[^and use it as the second criterion (it should be noted 
that this criterion is rarely relevant, since it is usually less 
constraining than the criterion of twice the mean particle 
separation). 

At each redshift bin, the respectively higher value 
among 2iimps and defined as: 


i?™ = max(i?Vif,2i?^p,) (7) 

is used as the minimum radius to obtain the abundance 
of voids. Thus, at each redshift we consider the most 
stringent of the two constraints, one taking into account 
the survey’s features (2i?rnps), the other depending on 
cosmology (ReuU~i)- 

In the next Section we present the Fisher formalism 
used in this analysis. 


IV. FISHER ANALYSIS 

The Fisher formalism is a method to determine a lower 
bound on the expected uncertainties with which cosmo¬ 
logical parameters can be determined. We rely on this 
formalism to estimate the power of surveys in constrain¬ 
ing cosmological models with void abundances. 

We discretize the survey into separate redshift bins and 
assume that the number of observed voids in any bin is 
a random draw from a Poisson probability distribution. 
That is, if the expected number of voids in the bin from 
the theory is A)., whose calculation was the subject of 
Section II, then the observed number of voids. No, has 
the probability distribution 


N^oe-N, 

p{No\No) = . (8) 

We assume that the drawings from two different red¬ 
shift bins are uncorrelated so that the probability to find 
No^i in the first bin, No ^2 in the second bin and so on, 
simply is: 


L{No^i, Noa, A^e,2; ••■) — 




a=l 


No,a\ 


(9) 


The Fisher information matrix is defined as: 


9 In L 9 In L 
dOm dOn 


( 10 ) 


where E stands for the expectation value and the 9i de¬ 
note the cosmological parameters. For the case of the 
Poisson distribution we will have: 


91nL _ / No^a dNo^a ,,,, 

~ ^ ^ 

Considering Equation using E [A^o,a] = A^e.a and 

the fact that for the Poisson distribution 


nNo,aNo,b] = E[No,a]Sab+E[No,a]E[No,b 
= Ne,adab + Ne,aNe,b, 


n2i 


it immediately follows [45] that: 


_ dNe^a dNe^a 1 
^ dOjn de„ No,a 

n. ’ 


(13) 


where the right hand side is calculated for a fiducial set of 
values of the parameters 0. We note that, when deviat¬ 
ing from fiducial cosmology, the inferred void radius will 
change. We consider the modification of the void radius 
when the adopted cosmology is not the fiducial one as: 
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FIG. 1: We represent the prediction from the Sheth and Van 
de Weygaert model matched to the simulation. The shaded 
areas show the effects of the variation of the value of (var¬ 
ied by 0.05 and 0.1 in the Figure) for the ACDM case. The 
marginalization over 5v allows to take into account the fact 
that the Sheth and Van de Weygaert fit is qualitative, and 
only roughly matches simulations due to the discreteness of 
the galaxy distribution as opposed to a smooth density field. 
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FIG. 2: Matching of the abundance of voids in simulations 
for models with different wq. The lines show the abundances 
from simulations for models with different wq (and fixed Wa)- 
The prediction from the Sheth and Van de Weygaert model 
are represented with dashed lines. 


V. SIMULATIONS 


.^Eul — .^Eul, true 


dUfid /(cLzdII) \ 

/(dzdn) J 


dUt 


(14) 


The variation of the radius of voids if a cosmology dif¬ 
ferent than the fiducial is assumed needs to be taken 
into account when performing the Fisher analysis, since 
it affects the derivatives of equation We note that 
ktrue and Vfid need to be in the same units. While cor¬ 
rectly taking into account this effect (considering equa¬ 
tion [ 14 ] when we calculate the numerical derivatives for 


Ne in equation 131, for the present analysis we are ne¬ 
glecting further effects (such as the Alcock-Paczyhski ef¬ 
fect), which should be quantified in future analysis with 
mocks based on Halo Occupation Distribution models 
(that promise to be more accurate when further details 
about Euclid and WFIRST will be available) and for the 
application of this method to data. 

The Cramer-Rao theorem states that any unbiased es¬ 
timator of the parameters 9a will have a covariance ma¬ 
trix larger than, or at best equal to, the inverse of the 
Fisher matrix. 

Thus the Fisher matrix formalism can give an esti¬ 
mate of the constraints to be obtained from cosmological 
probes in terms of area in the parameter space. 

While more sophisticated techniques exist to esti¬ 
mate such constraints [46], the Fisher matrix analy¬ 
sis is adapted for a first application of our simulation- 
calibrated number functions method. 


As discussed in Section I, we used simulations to tune 
the parameter 5y. We use simulations of sub-sampled 
dark matter and use the void finding technique described 
in |35] . While we leave for future work the tuning of Sy 
with mocks constructed using Halo Occupation Distri¬ 
bution models (which promise to be more accurate when 
further details about Euclid and WFIRST will be avail¬ 
able), the subsampling of simulations allows us to mimic 
the number density of tracers in the two surveys. 

For the simulations, we have used an adaptive treecode 
N-body method whose operation count scales as N log N 
in the number of particles, the 2HOT code |37|. With 
the use of a compensating smoothing kernel for small- 
scale force softening [48] and a technique to subtract the 
uniform background density, the error behavior and ac¬ 
curacy are particularly adapted for cosmological volumes 
(see [49] for details). To generate initial conditions we 
used a power spectrum calculated with CLASS m, and 
then realized the initial conditions with a modified ver¬ 
sion of 2LPTIC |ni| . The simulations assume the hducial 
Planck cosmology, we have then independently varied wq 
and Wa while adapting the amplitude of initial fluctua¬ 
tions such as to preserve the value of us determined by 
Planck. The size of the box is 1 Gpch~^ and it contains 
1024^ particles. We subsample the simulation to a mean 
density of h = 4 x 10“^ particles per cubic h~^Mpc, 
comparable to the number density expected from future 
surveys [13]. 

The comparison of the void abundance in the ACDM 
simulation with the theoretical model is shown in Figure 
[U The theoretical model is able to match the abundances 
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from simulations, except for smaller radius which we do 
not use to constrain parameters (as discussed in Section 
III). To parametrize the uncertainty in the determination 
of the parameter Sy we can marginalize over it to obtain 
the constraints. The shaded areas in Figure show the 
effects of the variation of the value of Sy (varied by 0.05 
and 0.1 in the Figure) for the ACDM case. The marginal¬ 
ization over Sy allows to take into account the fact that 
the Sheth and Van de Weygaert fit is qualitative, and 
only roughly matches simulations, due to the discrete¬ 
ness of the galaxy distribution as opposed to a smooth 
density field. 

Figures and show the comparison of void abun¬ 
dance with simulations for different values of the pa¬ 
rameters Wq and Wa- The used simulations mimic the 
number density of tracers in the two surveys, thus al¬ 
lowing us to match the abundance of voids in a realistic 
framework (although the volumes of the simulations are 
smaller than the actual volume to be obtained with both 
surveys). While at larger radii we are able to match the 
behaviour of simulations, we see that in Figures and 
[ 3 ] the theoretical model and the curves from simulations 
diverge at small radii. For voids of small radii, the void 
number function increases dramatically. As described in 
Section III, to ensure a robust estimation for the num¬ 
ber of voids, we cut off the constraints only considering 
voids at higher radii and thus avoiding to consider voids 
at smaller radii. 

Additionally, as mentioned before (and shown in Fig¬ 
ure]^, in order to account for the residual uncertainty in 
the Sheth and Van de Weygaert fit, we will marginalize 
over the parameter Sy. The abundance of voids obtained 
with the described methodology gives a measure of the 
power of the survey to constrain cosmological parame¬ 
ters, and is particularly robust, as based on recent ob¬ 
servations. The next Section presents the results of the 
calculation. 


VI. RESULTS 

The cosmological parameters we use to evaluate the 
Fisher matrix together with their fiducial values are given 
in table |l] (and correspond to Planck values [^). We 
use the BBKS analytical fit [53] for the transfer function 
T(k), and consider its generalization [^ that takes the 
effect of baryons into account. The matter power spec¬ 
trum at any given redshift can then be calculated as 


4 

Ps(k,a) = ^T^{k)D\a)-^-^Pcik) (15) 

Here, is the power spectrum of the primordial 
curvature perturbation and is assumed to be of the 

form A^(/c) = ^ 2 ^ 2 ^'^ = A^(fco)(A:/fco)"=*'^°^“^, where 
ko = 0.05 liMpc”^ is the pivot scale. Also, D{a) is 
the linear growth factor of structure, normalized so that 
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FIG. 3: Matching of the abundance in simulations for models 
with different Wa- The thick lines show the abundances from 
simulations for models with different Wa (and fixed wq). The 
prediction from the Sheth and Van de Weygaert model are 
represented with dashed lines. 


TABLE I: Parametrization of the cosmology and the hducial 
values chosen for the maximal set of parameters used in eval¬ 
uating the Fisher forecasts. 


fib 


h 

T 

flk 

Wo 

Wa 

Tis 

ln(10^°A(;(A:o)) 

0.049 

0.318 

0.67 

0.09 

0.0 

-1 

0 

0.96 

3.098 


at late time D{a = 1) = 1. In the case of a non-clustering 
dark energy fluid it can be found as the solution of the 
following differential equation: 


2 ^m(a) 

d In d 
”*"dlna 


d^ In d / d In d \ ^ 

d(lna)2 ydlnaj 

1 - ^ [^m(a) + {3w{a) + l)nDE(a)] 


(16) 


from [^ (to be more precise, D{a) is the growing mode 
of the solution of equation 16 we have S (x D{a)). 

To find out how a combination of different surveys can 
constrain the cosmological parameters, we simply need 
to add the Fisher matrices from those surveys. However, 
it is important to have the same set of cosmological pa¬ 
rameters for all surveys. We can change from one set 
of parameters to others by using F' = Jl.F.J in which 
Jij = dOildO'j is the Jacobian of the transformation. 

We consider Fisher forecasts computed from the 
Planck satellite (note that this forecast does not include 
CMB lensing); SN is the Fisher matrix for a Type la su¬ 
pernova survey such as the Large Synoptic Survey Tele¬ 
scope [171 [55] . 

Inverting the Fisher matrices, we can obtain the 1- 
a error ellipses in the wq — Wa plane, that measure the 
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capacity of a combination of experiments to constrain the 
dark energy equation of state. 

We show, respectively in Figures and the l-tr er¬ 
ror ellipses in the u>o — Wa plane for Euclid and WFIRST 
voids combined with Planck and HST (namely the Hub¬ 
ble Space Telescope Survey measurements of the Hubble 
constant from m)- For comparison, we also show the 
ellipse of Planck-|-HST alone, the Planck-l-HST-l- SN and 
the Planck -|- HST-I- Euclid BAO (from [H])- Figures]^ 
and thus allow us to compare the constraining power 
of void abundances for Euclid and WFIRST with BAO 
and supernovae. 

Additionally, to parametrize the uncertainty on the 
6y parameter, we can marginalize over it (as well as 
marginalizing over standard parameters). We vary the 
parameter 6^ by ±0.05, thus taking into account the ef¬ 
fects of sparsity (due to the discreteness of the galaxy 
distribution as opposed to a smooth density field) that 
might affect the void abundance and the ordering discrep¬ 
ancies between the simulations and the Sheth and Van 
de Weygaert fit seen in Figure For the marginalisa¬ 
tion to characterize the uncertainties and biases between 
the Sheth and Van de Weygaert model and the simula¬ 
tions we use a uniform prior (—0.45 ± 0.05). While we 
don’t expect our results to change substantially in the 
case of e.g. a gaussian prior (because the cosmological 
constraints are relatively robust with respect to the un¬ 
certainty in Sy), the precision in the calculation of the 
expected constraints is currently more strongly affected 
by the uncertainty about the details of the surveys than 
by the choice of the prior. Even though it is not the 
purpose of this paper to perform a full Markov Chain 
Monte Carlo exploration of the probability distribution, 
we point out that it would be the optimal way to explore 
the probability distribution function to be considered for 
future work, when details of the surveys will be updated 
and thus more reliable, allowing a sensitive improvement 
in the precision of constraints. The marginalization over 
6y also allows to take into account the residual effects of 
peculiar velocities on void abundance for voids of larger 
size, as well as further uncertainties related to the wa¬ 
tershed definition of voids in the tracer distribution m- 
Figures and also show (light-blue filled region) the 
l-cr error ellipses in the wq — Wa plane for the satel¬ 
lites after marginalization. We observe that the ellipse 
is only slightly thickened, showing that the cosmologi¬ 
cal constraints are relatively robust with respect to the 
uncertainty in Sy. 

Interestingly, we notice that the void abundance with 
WFIRST, while leading to similar constraints on the pa¬ 
rameter wq, constrains the Wa parameter slightly better 
than Euclid. Nevertheless, the conditional error on Wa 
from Euclid is lower - this means that if we combine 
constraints from void abundances with other probes that 
constrain effectively wq (such as supernovae), we can ex¬ 
pect the resulting combined constraints on Wa to be bet¬ 
ter for the probe combined with Euclid (as we will see in 



FIG. 4: 1-a error ellipses for several combinations of data sets. 
From outer ellipse to inner: Planck-l-HST, Planck+HST-|-SN, 
Planck-|-HST-|-Euclid Voids. The light-blue filled contour 
shows the result of the marginalization over the parameter Sy. 
For comparison, we also show in the striped-filled contour the 
BAO constraints from Euclid (also added to Planck-l-HST), 
from [T5| . 
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FIG. 5: l-cr error ellipses for several combination of data sets. 
From outer ellipse to inner: Planck-pHST, Planck+HST-|-SN, 
Planck-PHST-PWFIRST Voids. The light-blue filled contour 
shows the result of the marginalization over the parameter Sy. 
For comparison, we also show in the striped-filled contour the 
BAO constraints from Euclid (also added to Planck-pHST), 
from [T5| . 




Indeed, so far we have only compared the constraints 
from the void abundance to the contraints from super¬ 
novae, but additional constraints are obtained if super¬ 
novae are considered jointly with voids. Figure shows 
the comparison of WFIRST and Euclid void constraints 
when both are added to Planck -|- HST -|- SN. It should 
be noted that the axis range is smaller in this Figure. 

When independently added to supernovae, the void 
abundance from Euclid and WFIRST results in an en¬ 
hanced constraint, although, as expected from the previ¬ 
ous discussion, Euclid combined with supernovae is most 
constraining. Both WFIRST and LSST promise to bring 
new SN data. We can expect that, when added to the 
complementary probe from void abundance, these data 
increase the quality of constraints. While voids alone al¬ 
ready bring important information to constrain the dark 
energy equation of state, the different orientation of el¬ 
lipses - related to the different nature of probes as well 
as to the different systematics affecting the probes - for 
SN and voids will be a consistent asset to fully exploit 
the complementarity of the two probes. The resulting 
allowed wq — Wa parameter space is considerably reduced 
when voids are added to the supernovae probe, thus 
showing the power of void number counts to constrain 
the dark energy equation of state. 

TABLE II: Total number of voids forecasted for each survey. 


Survey 

Sky Fraction 

Total Number of voids 

Euclid 

0.36 

7.8x10® 

WFIRST 

0.05 

2.5x10® 


Also, it is worth notice that, while WFIRST will find 
less voids than Euclid (see Table [lT|), it still brings a com¬ 
petitive constraint as far as the void abundance is con¬ 
cerned, because it is sampling the cosmic web at smaller 
scales (higher density). 

The error ellipses of the two surveys cover different 
parts of the parameter space; thus — although it re¬ 
mains non trivial to combine cosmological constraints 
from the two data sets |SH1 — if jointly considered, they 
promise to yield improved constraints on the wq — Wa 
plane. The different constraints on the parameter space 
from the two surveys can be understood considering that 
they have galaxy number densities peaked at different 
redshifts (namely, the Euclid peak galaxy number den¬ 
sity will be at z ~ 0.75 — 0.85, while the WFIRST peak 
will be at z ~ 1.30 — 1.45). 

In addition to having different degeneracy directions, 
the Euclid and WFIRST surveys will also have different 
observational systematic errors due to their different ob¬ 
serving strategies. In general, a space-based slitless red- 
shift survey possesses depth variations across the field 
of view of the instrument due to varying image quality, 
spectral resolution, and coverage (chip gaps and defects), 
and over the area of the survey due to variations in sky 


brightness and stellar confusion. 

The WFIRST tiling strategy results in most galaxies 
being observed multiple times at widely separated posi¬ 
tions on the focal plane and at different zodiacal bright¬ 
ness levels, which will help to mitigate effects that de¬ 
pend on field position and sky brightness. However, its 
warm telescope results in a position-dependent thermal 
background that is highest near the center of the focal 
plane. 

Both surveys use multiple dispersion angles (imple¬ 
mented on Euclid by changing the grism and WFIRST 
by rotating the entire telescope). While we leave an in¬ 
vestigation of these effects on void statistics to future 
work, we note that comparison of the data sets in their 
region of overlap should help to constrain the systematics 
models for both surveys. 

As mentioned, because of correlations it would be non¬ 
trivial to combine Euclid and WFIRST constraints in 
the range where the two surveys overlap; nevertheless 
we can obtain a close estimate of the combined power 
of void abundances from the two surveys by considering 
that WFIRST is deeper that Euclid in the covered space, 
and Euclid is much wider. Thus as a simple estimation 
we can consider constraints from a survey as Euclid but 
with a smaller volume (i.e. subtracting the volume of 
WFIRST) and then adding the constraints of this survey 
to the constraints from WFIRST. We show the estimate 
of constraints obtained with such methodology in Figure 
1^ As expected, this estimate shows that combining voids 
from the two surveys would allow to enhance the quality 
of constraints. 

We showed that the large volume and deepness of the 
surveys allows to constrain effectively the wq — Wa param¬ 
eter space by using a phenomenological prediction for the 
number of voids with Euclid and WFIRST based on sim¬ 
ulations. To further dramatically improve the potential 
of the method, one could also consider voids at smaller 
radii which have not been used in our forecast - this 
would need a theoretical model able to predict the abun¬ 
dance more in detail also for smaller voids, and the use 
of simulations with larger volumes to confirm predictions 
[55] . Since the abundance of voids increases rapidly with 
decreasing radius, this could lead to a powerful increase 
in statistics and thus in even better constraints. The re¬ 
sults we obtain so far are able to show the potential of 
the method and serve as an initial probe of the power of 
these surveys in the framework of void abundance. 


VII. CONCLUSION 

In this paper we have considered the information 
gained by using the number of observed voids as a func¬ 
tion of redshift to constrain the equation of state of dark 
energy. The void number count is well performing in 
terms of improving constraints on w(z), in particular 
when combined with other probes. The advantage of 
void number counts is that the method does not need a 
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FIG. 6: Comparison of I-ct error ellipses for WFIRST and 
Euclid voids when the supernovae data is also considered 
(in each case the void abundance constraints are added to 
Planck+HST+SN data). We also show the constraints ob¬ 
tained by considering an estimate of the combination of voids 
from both surveys (see Section IV). (It should be noted that 
the axis range is smaller in this Figure.) 


separate expensive survey, since it can use data from the 
existing or planned redshift surveys to give extra infor¬ 
mation about the nature of dark energy. 

The method of observed cluster counts has been used 
to constrain the equation of state of dark energy (or, 
more generally, cosmological models, e.g. [51 IBDH53] ). It 
is worth noticing that the observed void counts promise 
to give complementary constraints to the cluster method. 
The first advantage of considering the void method is that 
the systematics will be different than in the cluster case. 
The estimation of the number of voids is based on the 
void’s observed “geometry” (radius), while the method 
for clusters relies on the use of the cluster “masses”. 

Clearly, there are different systematic errors involved 
in their estimation, and the confidence of the detection 
of a signal will be greatly enhanced when they are com¬ 
bined. Additionally, as mentioned, the measurement of 
voids does not need a different survey. To extract cosmo¬ 
logical information it is possible to directly use surveys 
already mapping the cosmic web for probes such as BAO 
— thus obtaining a complementary constraint. 

The second advantage is that the linear density thresh¬ 
old for void formation is different than the one for halos, 


which means we are considering a different probe. Fi¬ 
nally, the redshift dependences of the number of clusters 
and voids are different, since the mass scales probed are 
different. The abundance of clusters and voids are thus 
sensitive to the equation of state of dark energy in dif¬ 
ferent redshift ranges. For such reason the use of voids 
can give complementary and independent constraints on 
dark energy. 

We have compared two upcoming satellites: Euclid 
and WFIRST; and we have showed that their measured 
void abundance will constrain the parameter space 
wq — Wa differently, thus optimally sampling the cosmic 
web at a different scale and redshift. This results in 
joint constraints of an increased precision. Finally, the 
allowed wq — Wa parameter space is further reduced when 
void constraints are added to the supernovae probe, thus 
showing the power of the void number counts method to 
constrain the dark energy equation of state. 
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